The unitary three-body problem in a trap 
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We consider either 3 spinless bosons or 3 equal mass spin- 1/2 fermions, interacting via a short 
range potential of infinite scattering length and trapped in an isotropic harmonic potential. For a 
zero-range model, we obtain analytically the exact spectrum and eigenfunctions: for fermions all the 
states are universal; for bosons there is a coexistence of decoupled universal and efimovian states. 
All the universal states, even the bosonic ones, have a tiny 3-body loss rate. For a finite range 
model, we numerically find for bosons a coupling between zero angular momentum universal and 
efimovian states; the coupling is so weak that, for realistic values of the interaction range, these 
bosonic universal states remain long-lived and observable. 
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With a Feshbach resonance, it is now possible to pro- 
duce a stable quantum gas of fermionic atoms in the uni- 
tary limit, i.e. with an interaction of negligible range and 
scattering length a = 00 1]. The properties of this gas, 
including its superfluidity, are under active experimental 
investigation |2j. They have the remarkable feature of 
being universal, as was tested in particular for the zero 
temperature equation of state of the gas 0. In contrast, 
experiments with Bose gases at a Feshbach resonance suf- 
fer from high loss rates 0, 0, @ ' an d even the existence 
of a unitary Bose gas phase is a very open subject 0. 

In this context, fully understanding the few-body uni- 
tary problem is a crucial step. In free space, the unitary 
3-boson problem has an infinite number of weakly bound 
states, the so-called Efimov states In a trap, it has 
efimovian states 0,0] but also universal states whose en- 
ergy depends only on the trapping frequency . Several 
experimental groups are currently trapping a few parti- 
cles at a node of an optical lattice and are controlling 
the interaction strength via a Feshbach resonance. Re- 
sults have already been obtained for two particles per 
lattice node 0] , a case that was solved analytically . 
Anticipating experiments with 3 atoms per node, we de- 
rive in this Letter exact expressions for all universal and 
efimovian eigenstates of the 3-body problem for bosons 
(generalizing [tj to a non-zero angular momentum) and 
for equal mass fermions in a trap. We also show the long 
lifetime of the universal states and their observability in 
a real experiment, extending to universal states the nu- 
merical study of |10| . 

If the effective range and the true range of the interac- 
tion potential are negligible as compared to the de Broglie 
wavelength of the 3 particles, the interaction potential 
can be replaced by the Bethe-Peierls contact conditions 
on the wavefunction ip: it exists a function A such that 



^(ri,r 2 ,r 3 ) 
in the limit r?i = 
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taken for fixed positions 



of the other particle k and of the center of mass Ry of 
i and j. In the unitary limit considered in this paper, 
a = 00. When all the are non zero, the wavefunction 
tjj obeys the non-interacting Schrodinger equation 
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tp = Eip. 



(2) 



lo is the oscillation frequency and m the mass of an atom. 

To solve this problem, we extend the approach of Efi- 
mov H,Q| to the trapped case, and obtain the form 

^(n,r 2 , r 3 ) = ^ cm (C)F(R) (l + q) — V (a)Y{%p/p). 

(3) 

Since the center of mass is separable for a harmonic trap- 
ping, we have singled out the wavefunction -0 cm (C) of its 
stationary state of energy E cm , with C = (ri + r 2 + r3)/3. 
The operator Q ensures the correct exchange symmetry 
of ip: for spinless bosons, Q = P13 + P23, where Py trans- 
poses particles i and j; for spin 1/2 fermions, we assume 
a spin state TIT 80 that Q = — P 13 . The Jacobi coordi- 
nates are r = r 2 — ri and p — (2r 3 — I - ! — r 2 )/v3. Y" 1 
is a spherical harmonic, / being the total internal angu- 
lar momentum of the system. The function y(o:), where 
a = arctan(r//a), solves the eigenvalue problem 

- ip" {a) + V (a) = s 2 <p{a) (4) 
cos z a 

V (tt/2) = (5) 

v >(0) +r] (-l) l ^= ^(vr/3) = (6) 

with rj = —1 for fermions, 77 = 2 for bosons. An analytical 
expression can be obtained for J li ) . which leads to 
the transcendental equation for s [l6j: 
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FIG. 1: The constants s^ n for (a) 3 equal mass fermions and 
(b) 3 bosons, obtained by numerical solution of the transcen- 
dental equation Eq.(|7J. We have not represented the s; =0jn= o 
solution for bosons, which is purely imaginary. According to 
Eg. 1101 . each real Si t7l gives rise to a semi-infinite ladder of 
universal states. Note that the ground universal state has a 
total angular momentum I — 1 for fermions (E ~ 4.27/kj) 
and I = 2 for bosons (E ~ 5.32/kj). 



with the notation (x) n = x(x + 1) . . . (x + n). This equa- 
tion is readily solved numerically: for each I, the solu- 
tions form an infinite sequence (s;, n )n>0! see Fig^ As 
we show below, all solutions are real, except for bosons 
in the I = channel, where a single purely imaginary so- 
lution exists, s/ = o,n=o = so — i x 1.00624, the well known 
Efimov solution. Finally, the function F(R), where the 
hyperradius is R = J (r 2 + p 2 )/2, solves the problem: 

-|K^ + ^) + H F( * ) = (£ - £ ™ l)F<R) 

(8) 

where U(R) = h 2 s 2 /(2mR 2 ) + muj 2 R 2 /2, s being one 
of the si t n- This is the Schrodinger equation for a fic- 
tious particle of zero angular momentum moving in two 
dimensions in the potential U(R). 

When s 2 > 0, one takes s > and the solution is 



F(R)=R s e- R2 / 2 <L< q s \R 2 /al ) 



(9) 



where ctho = (h/muj) 1 ^ 2 is the harmonic oscillator length, 
Lq ' is the generalized Laguerre polynomial of degree q, 
q being an arbitrary non-negative integer. The resulting 
spectrum for the 3-body problem is 



E = E cm + (s Ln + 1 + 2q)hw. 



(10) 



The quantum number q leads to a semi-infinite ladder 
structure of the spectrum with a regular spacing 2ftw. 
This is related to the existence of a scaling solution for 
the trapped unitary gas ^3 an d the subsequent embed- 
ding of the Hamiltonian in a 50(2,1) algebra [3, 

lead- 
ing to an exact map ping between trapped and free space 
universal states |l9j. 



When s 2 < 0, as is the case in the I = n = chan- 
nel for bosons, the Schrodinger equation Eq.JSJ does not 
define by itself an Hcrmitian problem and has to be sup- 
plemented by a boundary condition for R — > |2fJ, |21| : 



F(R) oc Im 
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where R t is an additional 3-body parameter. For the 
resulting efimovian states, the function F is given by 

F(R) = R- 1 W [E _ Ecm)/2h ^ Sa/2 (R 2 /al ) (12) 

where W is a Whittaker function, and the energy solves: 

'l + s -{E-E cm )/huj 



argT 



-argT(l + sq) mod n. 



= -|s | ln(i? t /a ho ) 
(13) 



We did not yet obtain all the 3-body eigenstates |22j . 
Indeed, all the above states satisfy the contact condition 
(|T|) with a non-zero function A. But there are wavefunc- 
tions of the unitary gas which vanish when two particles 
are at the same point; these are also eigenstates of the 
non-interacting case. An example is the Laughlin state 
of the Fractional Quantum Hall Effect [2j| : 

i> = e-^ r ^ 2 < JJ [{x n + iy n ) - (x m + iy m )] M . 

1 < ri < m < 3 

(14) 

In the limit of high energies E 3> fkv, there are actually 
many of these A = states: their density of states (DOS) 
is almost as high as the DOS of the non-interacting case: 
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In contrast, the DOS of the A ^ states is only 
PA^o(E) 
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(15) 



(16) 



Ea. (|16|) is a consequence of Ea. (|17|) given be- 
low. We found Ea. l|15|l by applying the rank the- 
orem to the operator ^ (r 1 ,r 2 ,r 3 ) h-> (-0 o (ri, r l7 r 3 ), 
i>o( r i, r 2, ri), ipo(ri, r 2 , r 2 )) which associates, to each 
non-interacting eigenstate ipo of energy E, 3 functions 
of 2 atomic positions, and whose kernel is the space of 
A = states of energy E j^. 

This completes our derivation of all eigenstates of the 
unitary 3-body problem in a trap. Three types of states 
are obtained in general: universal eigenstates common 
to the non-interacting case, universal interacting states, 
and efimovian states depending on a 3-body parameter. 

We now prove that the Efimov effect is absent for 3 
equal mass fermions. This fact is known but to our 



knowledge not demonstrated. Numerically one can only 
check the absence of imaginary solution of the transcen- 
dental equation in some finite interval of s and I. Here 
we prove that for any I and any imaginary s, there is 
no solution to the problem (|4I5I6[) . Let us assume that 
s 2 < 1(1 + 1), and that I] 11511 are satisfied. We will show 
that the quantity Q(l, s 2 ) = ip'(0) - (-1) 1 (4/V3)<p(t./3) 
is non zero, which is incompatible with ©. We rewrite 
i fi} as f"(a) — u(a)ip(a). This is Newton's equation, 
a being the time and tp the position of a fictitious par- 
ticle subject to an expelling harmonic force with time 
dependent spring constant u(a;l, s 2 ) = ^t^j — s 2 > 0. 
Eq.© imposes that this particle reaches the origin at 
'time' 7r/2. The particle then should not reach the ori- 
gin earlier, otherwise the expelling force would prevent 
it from turning back to <p = 0. We thus can take the 
normalization <p(0) = 1, which implies <p'(0) < and 
ip(a) > for < a < tt/2. Thus, Q(l,s 2 ) < for I 
even. For I odd, one needs two intermediate results: (i) 
Q(l = 1, s 2 = 2) < (which we check by explicit calcula- 
tion); (ii) if (pi, f>2 are two solutions with u%>ux, then 
¥>2 < <Pi, and Q2 < Qi- because the spring constant for 
particle 2 is larger, particle 2 has to start faster and walk 
constantly ahead of particle 1 in the race towards the ori- 
gin, to satisfy Eq.JSJ). Now the assumption s 2 < 1(1 + 1) 
implies u(a;l,s 2 ) > u(a;l — l,s 2 — 2). One concludes 
that: Q(l,s 2 ) < Q(l = l,s 2 = 2) < 0. For bosons, we 
proved similarly that all the s 2 are positive, except for 
the well known s„ = o,;=o — i x 1.00624. 

It appears clearly on FigQ]that s; ; „ gets close to an 
integer value s^ n as soon as I or n increases, with 
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si, n = l + l + 2n for I > \r]\ 

si[ n = 2n-l + (2r) + ll)/3 for / < |r?| 



(17) 



To check this analytically, the transcendental equation is 
not useful. We rather applied semi-classical WKB tech- 
niques to the problem (|4I5I6[1 . and obtained |25|: 
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(18) 
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with /limax — 0.5357 the maximum of the Airy function. 

We now discuss the lifetime of the 3-body states found 
here in the trap, due to 3-body recombination to a deeply 
bound molecular state. The recombination rate is com- 
monly estimated as T\ oss oc Ph/(ma 2 ), where a is the 
range of the interaction potential, and P is the proba- 
bility that R < a 26] . Evaluating P from the 3-body 
wavefunctions obtained above for the zero range model, 
this gives for E not much larger than hio: 

2 s 
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loss 




(21) 



a/a ho 



FIG. 2: Numerical solution of the separable potential model: 
(a) 3-body eigenenergies and (b) predicted 3-body loss rates 
(for the case of 133 Cs, see text) , as a function of the potential 
range a (lower axis) and the 3-body parameter Rt (upper axis) 
(27) . (a) The lowest energy universal branch (stars) and an 
efimovian branch (crosses) have a very weak avoided crossing 
(inset). The analytical predictions of the zero-range model 
(solid lines) are in good agreement with the numerics (except 
for the avoided crossing); a linear extrapolation of the stars 
to c = matches the zero-range result at the 10 -3 level, (b) 
The universal states have a loss rate much smaller than ui. 



for a universal state with exponent s, and r^ s ™ oc u> 
for an efimovian state. Since s > 1.77 for fermions and 
s > 2.82 for bosons (FigQJ, Eq.J2U indicates that the 
lifetime of universal states is 3> l/u> for a <C aho- 

The existence of long-lived bosonic states is an unex- 
pected feature that we now investigate in a more realistic 
way. The unitary three-body problem in an isotropic har- 
monic trap may be realized experimentally by trapping 3 
atoms at a site of a deep optical lattice, and using a Fesh- 
bach resonance. For a broad Feshbach resonance, the ef- 
fective range is of the order of the Van der Waals length, 
which is roughly one order of magnitude smaller than 
aho, for an usual lattice spacing of ~ 0.5/im and a lattice 
depth of ~ 50 recoil energies. This experimental situa- 
tion is not deeply in the asymptotic regime of a zero range 
potential. Moreover, in the zero range model, there are 
energy crossings between universal and efimovian states 
as a function of Rt/a\ w (see solid lines in Fig|2t); as we 
shall see, for a finite range, there is a coupling between 
I = universal and efimovian states, leading to avoided 
crossings [i^, and to an additional contribution to the 
loss rate of I = universal states not included in Ea.l(2*Tjl. 

We therefore solve a finite interaction range model, the 
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Gaussian separable potential of range a |lfj , defined as 

(n.ralVKya) = - * 5 e^^/^ 2 5(R 12 -R' 12 ). 

(22) 

This leads to an integral equation that we solve numer- 
ically. In FigEK, we show two I = energy branches as 
a function of <j, corresponding in the zero-range model 
to the lowest I = universal state and to an efimovian 
branch. The smallness of the avoided crossing between 
the two branches shows that the coupling due to the finite 
range of the interaction is weak: the energy splitting at 
the avoided crossing is Ml ~ O.Olfkv, see inset of FigEt- 
We now revisit the calculation of the 3-body loss rate 
for bosons, since Ea. (|21|l neglects the contamination of 
the universal state by the efimovian state. To account 
for the losses we add to the Hamiltonian H scp of the 
separable potential model an antihermitian part leading 
to the effective Hamiltonian in second quantized form 

H cS = H scp -iB 3 ^ f [V f (r )] W)] 3 df (23) 
12m J 

where B3 is a numerical factor, whose actual value de- 
pends on short-range atomic and molecular physics. Spe- 
cializing to 133 Cs, we adjust the parameters of our model 
to B3 = 25 and a = 6.5nm in order to reproduce the 
three-body loss rate measured in a non-condensed gas 
for several negative values of a in |5j. To obtain the 
loss rates shown in Fig|2b, we restricted H e g to the two 
branches of Fig[2Ji: the eigenvalues of the resulting 2x2 
matrix have complex parts — ihTi oss /2. For the efimovian 
states, ri oss ~ 0.07w. For the universal states r loss is sev- 
eral orders of magnitude smaller; this remains true on the 
avoided crossing, because the coupling Q/2 of the univer- 
sal state to the efimovian state is much smaller than the 
decay rate of the efimovian state [29| . 

Experimentally, if one starts with the non-interacting 
ground state, a superposition of 3-body unitary eigen- 
states can be prepared by switching suddenly the scat- 
tering length from zero to infinity. The Bohr frequencies 
in the subsequent evolution of an observable would give 
information on the 3-body spectrum. For bosons, there 
will be a finite fraction of the sites where the three atoms 
have a long lifetime. This fraction is equal to the prob- 
ability of having populated a universal state, which we 
calculate to be ~ 0.174, a value dominated by the con- 
tribution (~ 0.105) of the lowest / = universal state. 

In summary, we obtained the complete analytical so- 
lution of a zero-range unitary 3-body problem in a trap. 
For bosons, there are efimovian and universal states, 
while for equal mass fermions we proved that all states 
are universal. All universal states are stable in the zero- 
range limit with respect to 3-body losses, not only for 
fermions, but also for bosons. From the numerical so- 
lution of a finite range model, we find that, although 
the bosonic universal states of zero angular momentum 



slightly mix with the efimovian states, their lifetime re- 
mains much larger than the oscillation period in the trap. 
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